Quantum Theory of a Resonant Photonic Crystal 
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We present a quantum model of two-level atoms localized in a 3D lattice, based on the Hopfield 
polariton theory. In addition to a polaritonic gap at the excitation energy, a photonic bandgap opens 
up at the Brillouin zone boundary. Upon tuning the lattice period or angle of incidence to match 
the photonic gap with the excitation energy, one obtains a combined polaritonic and photonic gap 
as a generalization of Rabi splitting. For typical experimental parameters, the size of the combined 
gap is on the order of 25 cm -1 , up to 10 s times the detuned gap size. The dispersion curve contains 
a branch supporting slow-light modes with vanishing probability density of atomic excitations. 

PACS numbers: 32.80.Pj, 42.50.Fx, 42.70.Qs 



I. INTRODUCTION 



Photonic crystals (PhCs)-systems in which the index 
of refraction varies periodically on the scale of light — are 
known to have an extraordinary ability to control the 
flow of lighljii^i^. While the underlying index of refrac- 
tion in conventional PhCs is often taken to be the bulk 
value, interesting effects can occur when the underlying 
medium possesses resonances at wavelengths comparable 
to the lattice spacing; for instance, one can dramatically 
widen the photonic bandgap by tuning the bandgap fre- 
quency (e.g. by changing the lattice period) to match the 
resonance frequency. Such "resonant PhCs" can be real- 
ized using cold atoms in optical lattices^iS, PhCs made 
from polaritonic materials^, and multiple-quantum-well 
array o 8 i 9 i 10 . Here, we will concentrate on the first class of 
resonant PhCs, originally analyzed by Deutsch, Spreeuw, 
Rolston, and Phillips^, who modelled the atoms in a ID 
optical lattice as a set of classical polarizable planes and 
showed that the interaction of the resonances with the pe- 
riodicity of the system gives rise to a photonic bandgap. 
Subsequently, van Coevorden et. al& extended this study 
to 3D by solving Maxwell's equations in a lattice of res- 
onating point dipoles using a t-matrix analysis. 

In this paper, we present a simple 3D quantum me- 
chanical model of an atomic PhC in which the elementary 
excitations are polaritons: coherent superpositions of 
atomic excitations and photons. Several features of pre- 
vious classical models appear naturally, and with some- 
what simpler interpretations, in the quantum model. For 
instance, we show that the resonance-induced bandgap 
arises as a generalization of Rabi splitting in a microcav- 
ity Our model also exhibits the important "Bragg reso- 
nant" modes first studied in ID by Deutsch et. ai, who 
identified them with the standing electromagnetic wave 
that supports the optical lattice^. Here, the Bragg reso- 
nant modes generalize to a family of modes occupying the 
boundary of the first Brillouin zone (BZ), and attached to 
the dispersion curve associated with atomic excitations; 
near the BZ boundary, they possess low group velocity 
but involve little excitation of the underlying medium, 
unlike states in "slow light" systems^. 



II. MODEL HAMILTONIAN 

Consider N localized two-level atoms of the same type 
in a fully-filled 3d cubic lattice, at sites r,; with lattice pe- 
riod t. To facilitate calculation, we enclose the lattice in 
a periodic electromagnetic cavity of volume V, which re- 
produces the physical behavior inside a sufficiently large 
lattice. The Coulomb-gauge Hamiltonian is 
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where e is the energy difference between the atomic levels, 
b\ = |l)(0| i and bi = |0)(l| i are the level raising and low- 
ering operators for atom i, and at and ag ff are creation 
and annihilation operators for photons with wavevector 
k and polarization a . A(r) is the vector potential, 
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where e^ a is the unit polarization vector for at . 

Let us suppose that the average number of atomic ex- 
citations in the system at any time is much less than N. 
In that case, the atomic excitations are approximately 
bosonic, in the same sense that spin waves are bosons^ 3 -. 
Therefore, the two photon polarizations, which excite or- 
thogonal atomic states, decouple for each k. We thus 
drop the a label, with the understanding that the dis- 
persion relations we will later obtain are doubly degen- 
erate. This is also consistent with the weak polarization 
dependence obtained by van Coevorden et. al& In con- 
trast, polarization effects play an important role in con- 
ventional PhCs 3 , as well as multiple-quantum- well reso- 
nant PhCs&i^, due to the finite size of the scattering 
centers. 

We can treat the in (Q]) as numbers (perfect lat- 
tice positions) rather than operators, since the electronic 
wavefunctions are typically much narrower than the lat- 
tice spacing. At each site, the momentum operator is 

Pi = - mexoiipl - bi) , x QX = (l|sf)0). (3) 
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Let us also define momentum-space excitation operators 
-i=£^^H4fE^ (4) 



where the wavevectors q are restricted to the first BZ, 
corresponding to the fact that an excitation "wave" has 
no meaning between lattice points. As explained above, 
these are approximately bosonic: [bg, bt,] ~ 5^. 

Substituting ©-jH) into ([1]), we obtain the microscopic 
polariton Hamiltonian first derived by HopfielcU^ in the 
context of crystalline solids: 
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where G's run over all reciprocal lattice vectors, and 
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where a is the fine structure constant. The atom-photon 
interaction consists of two parts. The first part, on the 
second line of (0, describes the lattice absorbing a pho- 
ton with wavevector q + G to create an atomic excitation 
with wavevector q, and the reverse process of destroying 
an excitation to emit a photon. The remaining inter- 
action terms describe the creation and annihilation of 
associated pairs of photons and atomic excitations. The 
usual way to diagonalize ([5]) is to introduce polariton 
operators a ^ 13 ' 14 for each reduced wavevector q, as lin- 
ear combinations of b^, a^ + g, and a_^^ (for all 
g). Stipulating that these act as decoupled lowering op- 
erators for H, one obtains the polariton energies as so- 
lutions of a (2n + 1) x (2n + 1) eigenvalue problem for 
each (f, where n is the number of BZs included in the 
calculation. Higher BZs were first included into the Hop- 
field theory by Knoester and Mukamel~ in their calcula- 
tion of polariton-mediated intermolecular forces in solids. 
There, the photons in the higher BZs were taken to be 
decoupled from the atomic excitations, which was appro- 
priate since the BZ energy was many orders of magnitude 
larger than e. In our system, the two energies are com- 
parable, and we must incorporate the interaction up to 
at least the second-order zones. 

It simplifies the calculations to drop the "counter- 
rotating" interaction terms in ([5]) describing the creation 
and annihilation of pairs. This is physically justifiable 
even though the discarded terms have the same cou- 
pling strength C^ + g as the remaining interaction terms, 
because the pair creation and annihilation process is a 
quantum mechanical fluctuation of the "vacuum" with a 
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FIG. 1: (Color online) Single-polariton dispersion for a 3d 
cubic lattice along [100] in the extended zone scheme, with 
e = 3eV and three different coupling strengths, associated 
with the parameters xai = (non- interacting), lA(C Q ~ 
0.18 meV), and 2 A (C Q ~ 0.35 meV). The vertical dashed line 
indicates the BZ boundary at |q| = 1.00025 e/hc. The graphs 
are generated numerically from 0, summing over 125 BZs. 



finite energy gap e + 7ic|g|. For e and hc\q\ both on the 
order of eV, and lattice periods at optical wavelengths, 
Cq ~ 10 -4 eV <C e + Tic\q\. Such fluctuations are thus 
extremely rare and have a negligible effect on particle 
energies. The interaction terms describing the conver- 
sion of a real photon into an atomic excitation, and vice 
versa, remain important: since the existing particle pos- 
sesses energy, these processes involve a much smaller en- 
ergy fluctuation. The approximation holds provided we 
look at values of \q\ comparable to both e/hc and the BZ 
energy, which is exactly the regime we are interested in. 

The Hamiltonian now decouples into N independent 
pieces, H = ^2gHg, one for each reduced wavevector: 



Hq = eblb q ~+J2 hc W+G\ 4+5 a <f+S 

G 



(7) 



This says that each photon mixes with all other photons 
having wavevectors that differ by a reciprocal lattice vec- 
tor, as one expects of a PhC system. Here, the mixing is 
mediated by the atom-photon interaction. Since |7| has 
the quadratic form . PjH-ij (3j , it can be diagonalized 
as J2 n EnCtti a n, where the cv's are boson operators de- 
fined by a n — J2j w j n flj> En is the nth eigenvalue of H, 
and w n is the corresponding eigenvector. We can thus 
obtain the polariton energies E'S by including a finite 
number of BZs in the sum and diagonalizing the associ- 
ated matrix. 



III. BAND STRUCTURE 

Fig- [H shows the polariton dispersion curves along the 
[100] direction for a blue-detuned optical lattice. The in- 
teraction opens up two energy gaps in the polariton spec- 
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trum: an indirect "polaritonic gap" A po i at e due to the 
repulsion between the bare dispersion curves, and a pho- 
tonic bandgap A p b g at 7ic|q| where Q is the BZ boundary. 
We have also calculated the density of polariton states; 
after integrating over all angles, we find that the density 
of states is enhanced near the band edges, but remains 
nonzero at all energies because the exact sizes and posi- 
tions of the gaps vary with angle. The system therefore 
does not possess a complete gap, essentially because of 
the weakness of the electromagnetic interaction. The gap 
sizes vary continuously as we change the lattice period a, 
and thus V (keeping N and all other parameters con- 
stant). As shown in Fig. [2 the gaps meet and become 
significantly enhanced when the BZ boundary intersects 
the crossing point of the bare dispersion curves. 

To understand the nature of the spectrum at the BZ 
boundary, consider a photon with wavevector k = Q 
along one of the faces of the cube. There is another such 
photon, with wavevector Q + G lying on the opposite 
face, such that |q| = |q + G |. (When Q lies on an edge or 
corner of the BZ boundary, there are more partners; we 
will not consider these cases, but they can be treated in 
a similar fashion.) The two photons mix strongly since 
they have the same energy, so we can neglect the other 
photon states and use the effective Hamiltonian 

t r 
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Thus, the polariton energies at the BZ boundary are 
hc\Q\, 

hc\< 



± 



Hc\q\ 



2Ct 



(9) 



These are exactly the energy levels resulting from Rabi 
splitting of a two-level atom interacting with two coun- 
terpropagating photon states with wavevectors ±Q, with 
an effective cavity size V/N. In the exactly-tuned case 
e = 7ic|q|, has a special significance: as shown in 
Fig. Hfb), these are the upper and lower edges of the 
bandgap. The resonant enhancement of the bandgap in 
this system is thus a manifestation of the Purcell effect 1 -^. 
Intuitively, we can imagine enclosing a single atom in a 
microcavity with the dimensions of the unit cell; if the 
cavity walls are mirrors, the atom sees a lattice of atoms 
similar to the one considered here. 

We have checked © against numerical solutions of ([7]) 
including the 125 lowest BZs, for various values of Q along 
the BZ boundary up to 40° from the [100] direction. For 
e = 3eV and xoi = 2 A, the error is always less than 
0.02 cm -1 , three orders of magnitude smaller than the 
maximum gap size. 

The size of the gaps in the exactly-tuned limit can be 
estimated by substituting e — hc\Q \ into ([9]): 
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FIG. 2: (Color online) Single-polariton dispersion along [100], 
with e = 3eV, xoi = 2 A, and different lattice periods: (a) 
|q| = 1.00025 e/hc, and (b) |q| = e/hc. Plots (c) and (d) 
show the corresponding overlaps of the polariton with the bare 
excitation, (0| 6^at|0), for the polaritons on the dispersion 
curve leading to the purely photonic state at q — Q (indicated 
with arrows in (a) and (b)), which have no atomic component. 



For e, 7ic|Q| ~ 3eV, and x 01 fa 2 A, A/hc w 25 cm -1 
(~ 10 -4 e), in agreement with Fig. [2jb). We can also ob- 
tain limiting expressions for the gaps when they are sig- 
nificantly decoupled. Consider |q| > e/hc, as in Fig. [^a). 
Away from the BZ boundary, we can neglect the effect 
of photons in higher BZs, and the effective Hamiltonian 
matrix is H = [e, —iC^; iC^, 7ie|<f|], with eigenvalues 



± = e + fic|g1 



hc\q\ 



(11) 



The contribution to the indirect polaritonic gap from the 
large-g branch of the dispersion curve, which is truncated 
at the BZ boundary, is obtained from the large-q expan- 
sion of Ijlip evaluated at q = Q. The contribution from 
the small-g branch cannot be found by setting q = in 
(|11|) due to our preceding approximations, so we instead 
calculate an upper bound on it by evaluating it at the 
minimum, \q\ = e/2hc. The resulting polaritonic gap is 
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With the same lattice parameters, A pol sa 10 3 cm 1 (~ 
10 -8 e). From the large- |q| expansion of ©, the photonic 
bandgap is A pbg = C%/%c\q\, strictly smaller than (fT2"j) . 
Therefore, the effects of the polaritonic interaction are 
very small when the system is detuned. 

This model can also be used to study the quasi-lD ge- 
ometry considered by many authors, in which atoms are 



trapped along periodically-stacked infinite sheets. Con- 
sider a 3d lattice in which the lattice spacing in one of 
the directions, i\, is much larger than the spacing in the 
other two directions. The relevant wavevectors, lying on 
the BZ boundaries closest to the origin, have magnitude 
|Qx| = 7r/^i and point in the direction of stacking. In this 
regime, this model can be directly compared with the 
semiclassical analysis of Deutsch et. al. 5 . For instance, 
the semiclassical theory predicts bandgaps from E [d) to 

e and from ?ic| Qi | to E^ 1 ' for blue-detuned lattices. A 
short calculation, using Eq. 15-19 of that paper, yields 

{d) _ e + hc\Q\ \( e — ^.c|Q| \ 2 3h 2 cyrj 

± ~ 2 V V 2 J ' 2|Qi| 1 ' 

where r] is the surface density along each sheet and 
7 <C (E± — e)/h is the linewidth of the atomic transi- 
tion. Using the golden rule prescription for the natu- 
ral linewidthii, 7 = (Acte 3 x^) / (3h 3 c 2 ) , this reduces to 
© with C 2 replaced by C 2 • e/hc\Qi\. The bandgaps 
predicted by the semiclassical and quantum mechanical 
theories are thus similar for e ~ Tic | Q 1 1 , which is also 
the regime where the bandgaps are significant. In the 
exactly-tuned case, the results are identical, and one ob- 
tains 

A ld = 2ex 01y ^j. (14) 

Actual Id/ 2d lattices are more problematic since each 
atomic excitation is coupled to photons with a continuum 
of wavenumbers in the transverse direction, which smears 
out the gaps. One might avoid this using an actual cavity 
in the transverse direction, making the electromagnetic 
field effectively 1d/2d. 

IV. SLOW POLARITON MODES 

The energy E® in © corresponds to a polariton cre- 
ated by the operator (<Xq — ol +G r )/"\/2- This remains an 
exact polariton state when we include higher BZs in the 
effective Hamiltonian. (In fact, there is a family of such 
states for each pair of BZ boundaries.) These "purely 
photonic" polaritons are reminiscent of "dark states" in 
electromagnetically induced transparency (EIT)i£, since 
§8§ is identical to the EIT effective Hamiltonian with the 
atomic excitation and two photon modes acting as the 
levels of the A system. In EIT, a "dark state" arises: 
a coherent superposition of atomic levels that does not 
couple to the radiation. The analog in our case is a 
non-interacting photonic state, with no atomic compo- 
nent. Its classical limit is a standing electromagnetic 
wave commensurate with the lattice. Since the laser light 
that supports the lattice always falls exactly on the BZ 
boundary^, the stability of the optical lattice relies on 
the existence of such standing wave modes; other modes 
are Bragg reflected away. In a sense, the lattice "se- 
lects" the standing wave modes from the incoming laser 
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FIG. 3: (Color online) Photonic gap at wavevectors Q = 
[k x ,7r/£, 0] along the BZ boundary, for e = 3eV, X01 = 2 A, 
and tt/£ = 0.9 e/hc = 1.4 x 10 5 cm -1 (red-detuned). The 
dashed lines show k y vs. k x for the surface \k\ — e/hc and 
the BZ boundary; here, the ordinate is not drawn to scale. 
The gap is largest at the intersection of the two surfaces, 
i.e. |q| = e/hc. 

light. Similar modes have been observed in other res- 
onant PhC system o 8 ' 9 ' 10 . We have shown here that in 
the self-consistent limit of complete quantum coherence 
and low excitation density, this selection takes place at 
the quantum state level. Only the purely photonic po- 
laritons can support a macroscopic population, since they 
are the only elementary excitations of the interacting sys- 
tem with zero atomic component. 

In the 3D system, there is a family of purely pho- 
tonic polaritons everywhere on the boundary of the first 
BZ. Remarkably, these states are attached to the slow, 
"atomic" branch of the dispersion curve. These appear to 
be analogs of the slow, non-degenerate, longitudinal elec- 
tromagnetic modes that appear in the classical t-matrix 
calculation of Coevorden et. al. 6 Our model shows that 
the photonic component along this branch goes continu- 
ously from nearly zero to unity as we approach the BZ 
boundary, as shown in Fig. [^c) and (d). Therefore, by 
exciting polaritons over a range 

^^-^-lOcm- 1 (15) 
he he 

around wavevector Q, one could create a wavepacket that 
propagates slowly but has low atomic excitation density. 

V. CONCLUSION 

We have presented a quantum model for an atomic 
lattice that applies directly to optical lattices filled with 
cold atoms, containing behavior similar to other resonant 
PhC systems. The system possesses two gaps (polari- 
tonic and photonic) at each angle, and can be tuned so 
that the gaps meet to create a combined gap orders of 
magnitude larger than the individual detuned gaps, in 
a process analogous to microcavity Rabi splitting; how- 
ever, there does not exist a complete gap. The quantum 
analysis yields a branch of the dispersion curve that has 
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low group velocity and atomic component vanishing at 
the BZ boundary. 

These effects could be explored with alkali atoms held 
in a cubic lattice made by near-IR light, by introduc- 
ing a probe beam at an angle to the axis of the lat- 
tice. One should choose an atomic transition e such 
that 1 < ei/nhc < y/3, where I is the lattice period, 
and use probe wavevectors with magnitude lying in a 
range A/hc ~ 10cm _1 around |g| = e/hc, at an an- 
gle cos -1 (nhc/ei) to a lattice axis (Fig. [3]). Although 
the present theory applies to an infinite lattice, the pre- 
dicted frequency shifts may be observable close to the 
atomic resonance, even in a lattice of about 100 atoms 
on a side. 

We have treated the atomic positions as fixed, as would 
be the case for a strongly-confining optical lattice where 
the rate at which each atom tunnels to a different lattice 
site is negligible compared to the radiative lifetime. The 



presence of non-zero hopping amplitudes would add an 
imaginary part to the polariton energies, proportional to 
the tunneling rate. The size of the band gaps would be 
reduced by the corresponding amount. 

Finally, it is interesting to note that the gap in Eq. 
(|10p. which scales as e relative to the photon energy, is 
O(10 -2 ) e for X-rays. Aspects of this theory might thus 
be applicable to crystalline solids in the X-ray regime, 
where a similar effect — superradiant scattering enhance- 
ment due to nuclear resonances — is known to exists. 
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